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ABSTRACT 

Sizes of galaxies are an important diagnostic for galaxy formation models. In this study I use the abundance 
matching ansatz, which has proven to be successful in reproducing galaxy clustering and other statistics, to 
derive estimates of the virial radius, ^2oo> f° r galaxies of different morphological types and wide range of 
stellar mass. I show that over eight of orders of magnitude in stellar mass galaxies of all morphological types 
follow an approximately linear relation between half-mass radius of their stellar distribution, ri/2 and virial 
radius, r\ji as 0.015/?200 with a scatter of m 0.2 dex. Such scaling is in remarkable agreement with expectation 
of models which assume that galaxy sizes are controlled by halo angular momentum, which implies r\ji oc 
AR200, where A is the spin of galaxy parent halo. The scatter about the relation is comparable with the scatter 
expected from the distribution of A and normalization of the relation agrees with that predicted by the model 
of Mo, Mao & White (1998), if galaxy sizes were set on average at z ~ 1-2. Moreover, I show that when 
stellar and gas surface density profiles of galaxies of different morphological types are rescaled using radius 
r n = 0.015^200' the rescaled surface density profiles follow approximately universal exponential (for late types) 
and de Vaucouleurs (for early types) profiles with scatter of only as 30 - 50% at R » 1 - 3r„. Remarkably, both 
late and early type galaxies have similar mean stellar surface density profiles at R > lr„. The main difference 
between their stellar distributions is thus at R < r„. The results of this study imply that galaxy sizes and radial 
distribution of baryons are shaped primarily by properties of their parent halo and that sizes of both late type 
disks and early type spheroids are controlled by halo angular momentum. 



1. INTRODUCTION 

In the standard hierarchical structure formation scenario 
galaxies form at the minima of potential wells formed by non- 
linear collapse of peaks in the initial density field. The pro- 
cess of galaxy formation is expected to be complex, highly 
nonlinear process, involving forces on a wide range of scales, 
supersonic, highly compressible and turbulent flows of gas, 
and a variety of cooling, heating, and feedback processes. 

Despite the apparent complexity of formation processes, 
observed galaxies exhibit a number of tight scaling relations 
between their structu ral parameters and, as first shown b y Fall 
& Efs tathioul ( |1980J> and elaborated by Mo, Ma o, & White 
1998 1 see~a l so IDalcanton et aL] (|1997|; |Av ila-Rees e et al. 
1998 20081; IDutton et al.| H2007) ; |Dutton| ( |2009D ; |Fu et al. 
(2010)), these scaling relations can be reproduced in a fairly 
simple framework, in which sizes of galaxies are determined 
by the sizes of their initial rotationally supported gaseous 
disks, which, in turn are set by the angular momentum of 
gas. Under assumption that gas angular momentum is propor- 
tional to that of dark matter, such models predict that galaxy 
size should scale as oc AR2<m, where /?20o is the "virial" ra- 
dius defined as the radius enclosing overdensity of 200 with 
respect to the critical density of the universe, p cr (z), so that 
M 200 = (4;r/3)200p cr (2)4 )0 . 

In addition, regularity in galaxy properties is implied by 
success of the abundance matching ansatz, in which relation 
between total mass of halos, M, and stellar mass of galax- 
ies they host, M, is derived from a simple assumption that 
the relation is approximately monotonic and cumulative abun- 
dance of galaxies with masses above a given M, is matched 
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to cumulative number density of halos with masses above M. 
This model is remarkably successful in reproducing cluster- 
ing of galaxies of different l uminosities and at different red- 
shifts ( Kravtsov Tt al.||2004| |Tasitsiomi et al.||2004| |Conroy 
l et al.|2006| |Reddick et al.|2012|> and other statistics ( Vale & 



Ostriker 2004, 2006, Behroozi et al. 2010 



2010| |Moster et al.|2dl2||Hearin et al.|2012) . 



2012| iGuoetaL 



In this paper, I use the abundance matching ansatz to exam- 
ine relation between sizes of stellar systems of galaxies, char- 
acterized by the three-dimensional half-mass radius, r\i% and 
virial radius of their halos, Z? 2 oo> derived using the abundance 
matching ansatz. I show that over the entire observed range of 
stellar masses and morphologies, galaxies exhibit an approxi- 
mately linear scaling relation between stellar half-mass radius 
and halo virial radius with normalization and scatter consis- 
tent with expectation of the |Mo et al.| ( 1998) model. Further- 
more, I show that stellar and gas surface density profiles of 
galaxies rescaled using radius r„ = 0.015/?200 follow univer- 
sal profiles with a scatter as low as as 30-50% at intermediate 
radii within optical extent of galaxies. 

Throughout this paper I assume a flat ACDM model with 
parameters On = 1 - O a = 0.27, Q b = 0.0469, h = 
i7o/(100 km s~' Mpc" 1 ) = 0.7, cr 8 = 0.82 and n s = 0.95 com- 
patible with combined c onstraints from WM AR BAO, SNe, 
and cluster abundance (Komatsu et al. 2 01 

2. ABUNDANCE MATCHING 

To estimate the virial masses and radii of halos hosting 
galaxies, I use the abundance matching ansatz, in which rela- 
tion between total halo mass, M, and stellar mass of galaxies 
they host, M», is established implicitly by matching cumula- 
tive stellar and halo mass functions: n^(> M) = n g (> M„). 

A number of estimates of the M* - M relation using this 
technique has been presented in the re cent literature (e.g., 
|Moster et aT|2012| |Behroozi et al.|2012) . However, the re- 
lations derived in these studies are based on stellar mass func- 
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tions (SM Fs) known to under estimate abundance of massive 
galaxies (Bern ard! et al.|2010|l and the double power law fit to 
the M* — M relation of Moste r et al. ( |2012| l does not capture 
the upturn in the relation at M* < 10 M Q originating fr om the 
steepening of the stellar mass function at t hese masses (Baldry 
et al |2008[|2Tri2{IPapastergis et al.|20l2) . Therefore, in this 



study I re-d erive the M. - M rel ation to fix these problems 

I use the Tin ker et al.| ( |2008| l calibration of the halo mass 
function for M200, which was calibrated using host halos only. 
To account for subhalos, I correct the host mass function by 
subhalo fraction, / sub (> M) = [w tot (> M)-n host (> M)]/n host (> 
M), to get n tot (> M) - the mass function that includes both 
hosts and subhalos. The latter was calculated using current 
M200 masses for hosts and corresponding masses at the ac- 
cretion ep och for subhalo s using z — halo catalog halo 
catalog of Behroozi et al. (2011) derived from the Bolshoi 
simulation (Klypin et al. 20ll|) of (250/z 'Mpc) 3 volume in 



the concordance cosmology adopted in this study. The sub- 
halo fraction in the Bolshoi simulation is parametrized as 
/sub = min[0.35, 0.085(15 -log 10 M200)]. The halo mass func- 
tion derive d from th e Bolshoi simulation agrees within 5% 
with the |Tinker et al.| ( [2008[ ) parameterization, but the latter is 

more accurate at the highest halo masses. 

I combine tw o re cent calibrations of th e SMF by Pa paster-] 
gis et al. ( 2012) > and |Bernardi et al. ( 2010[ ) to accurately char- 



actenze SMF behavior at both small and large M*, respec 
tively. I use these two calibrations to construct a combined 
stellar mass function, n(M„) - max[«pi2, «bioL that spans 
from M* « 1O 7 M to M* « 1 12 M G . Both stellar mass func- 
tions assume Chabrier (2003) IMF to estimate stellar masses 
of galaxies. For ripn I adopt double Schechter form given 



by eq. 6 of Baldry et al. (2012 1 with the following param- 
eters: log 10 M. = 10.66, 0^ 3.96 x 10~ 3 , oti = -0.35, 
<p* 2 = 6.9 x 10~ 4 , 0L2 = -1.57. These parameters are in general 
agreement with the bes t fit parameters deri ved for the local 
stellar mass function by Baldry et al. ( 2012| ). Note that SMF 
at M* < 10 8 M Q is quite uncertain due to incompleteness of 



low s urface brightness galaxies in this regime ( |Baldry et al. 
|2012| ); the current SMF measurements at these stellar masses 
should be considered as lower limits and the actual SMF may 
be somewhat steeper still. For «bio I use parameter values 



given in the bottom row of Table 4 in Bernardi et al. (2010 
unbracketed values) and the Schechter parametnzation of the 
SMF given by eq. 9 in that paper. I refer readers to the origi- 
nal papers for further details on how the stellar mass functions 
were estimated. 

3. GALAXY SAMPLES 

To estimate the size-virial radius relation, I have selected 
several publicly available datasets chosen to span the entire 
range of galaxy stellar masses^] and morphologies. First, 
I use a compilation of stellar masses and effective radii 
for sp heroidal, early-type galaxies from Misgeld & Hilker 
(|201 \\. These include ellipticals (Es) and dwarf elliptical (dE) 



galaxies in the Virgo cluster with HST (Ferrarese et al.|2006]) 
and the VLT/FORS1 obs ervations of dEs in t he Hydra I and 
Centaurus clusters (Misg eld et al.|2008|[2009| ), and the dwarf 
spheroidal (dSph) galaxies in the Local Group. The sample of 
late ty pe galaxies includes the THINGS/HERACLES galax- 
ies of |Leroy et al.] ( [2008] ) and the LITTLE THINGS sa m- 
ple of dwarf irregular galaxies from Zhang et al. ( 2012[ ). I 



4 Stellar masses in all of the samples were estimated assuming the |Chabrier| 
IMF. 
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Fig. 1 . — Relation between half-mass radius of stellar distribution in galax- 
ies of different stellar masses (spanning more than eight orders of magnitude 
in stellar mass) and morphological types and inferred virial radius of their 
parent halos, K200, defined as the radius enclosing overdensity of 200p cr , and 
estimated as described in §12] The red pentagons and hexagons sho w a sam- 
ple of elliptical and dwarf elliptical galaxies from the compilation of Misgeld 
& Hilker 1 201 1 1; blue circles are the late typ e galaxies from the samples ot 
|Leroy et al.ltZoOT 



rraosia 

as described in § IS] while the star symbol shows the Milky Way; the red 
cartwheel poin ts show the Local Group d warf spheroidal galaxies from the 
compilation of Misgeld & Hilker (20111. The light blue dashed line and 
dot-dashed orange line show the average relations derived for late a nd early- 
type g alaxies, respectively, from the average S1/2 - M, relations of |Bernardi| 
|et al.H 2012t. Dark red shaded band shows 2<x scatter around the mean re- 



lation calculated for all galaxies in the |Bernardi et al.| (20 1 2 1 sample. The 
orange dot-dashed line with error bars shows the mean relation and 2<x scat- 
ter for massive SDSS galaxies presented in Szom oru et al.|(2012| ; individual 
galaxies from this sample are shown by blue (Sersic index n < 2.5) and red 
(n > 2.5) dots. The gray dashed line shows linear relation n/2 = 0.015i?200 
and dotted lines are linear relations offset by 0.5 dex, which approximately 
corresponds to the scatter in galaxy sizes from distribution of halo spin pa- 
rameter A under assumption that ri/2 °c ARiqq. 

also include the stellar mass profile of the Milky Way us- 
ing a combination of thin an d thick s tellar disks with pa- 
rameters given in Table 2 of McMillan (2011). For the late 
type samples, I used the deprojected stellar surface density 
profiles presented in these studies to estimate the half mass 
radius, r\ji, directly from profiles. The radius r\ji was de- 
termined as the radius that contains half of the stellar mass 
of galaxies using the cumulative mass profile of each disk: 

MJ<R) = 2n £l,(R')R'dR'. 

In addition, I use the average relations between half-light 
radius and stellar mass, (ri/2|M,)> derived for earl y and late 



type galaxies in the SDSS from the recent study by Bernardi 



|et al.| ( |2QT2| SerExp values in their Table 4). I also use intrin- 
sic scatter about the mean relation calculated for both early 
and late type galaxies (M. Bernardi 2012, priv. communica- 
tion). Finally, I use a half-mass radii and stellar masses for a 



sample of massive SDSS galaxies presented in Szomoru et al. 
pOlIl see their Table 1). 
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4. RESULTS 
4.1. The size-virial radius relation 

To derive the size-virial radius relation, I first assign M200 
to galaxies using their stellar mass and the M„ — M200 re- 
lation derived using abundance matching. I then estimate 
the three-dimensional half-mass radius from the projected 2D 
half-mass radii reported for observed galaxies. 

I assume that in late type galaxies stars are in a disk and 
hence the 2D R\^ radius is equal to the 3D r\/2 radius. For 
early type galaxies I assume that stars have spheroidal dis- 
tribution and convert projected R s into the 3D half-mass ra- 
dius using r\/2 = l.34R e . This expression is accurate for 
spheroidal systems described by the Seisic profile with a wide 



10 



range of the Sersic index values (see eq. 21 in Lima Neto et al. 
[19991 

Figure |T| shows the derived relation between ri/2 and the 
virial radius of parent halos, /?200> for all the observational 
samples described above. Remarkably, ri/2 scales approxi- 
mately linearly over two orders of magnitudes in radius and 
over eight orders of magnitude in stellar mass from the dwarf 
spheroidal galaxies to the most massive ellipticals, as ex- 
pected in the model s assuming that size of galaxies is pro - 
portional to AR 20 o ( |Fall & Efstofluou|l980l |Mo et al.p998) . 
The linear relation, r\/2 = 0.015/?200 is shown by the gray 
dashed line. The normalization of this relation is chosen so 
that distribution of points is approximately symmetric around 
the line. The formal power law fit to the ri/2 - Rim relation 
of individual galaxy points shown in the figure gives slope 
of 0.95 ± 0.065, normalization of 0.015 ± 0.0007, and scatter 
of 0.20 ± 0.016 dex (all errors indicating 95% confidence in- 
terval). These values should be taken only as indicative that 
relation indeed has the slope close to linear, given that galaxy 
sample only samples the mass range in a semi-random fash- 
ion. More robust determination of the slope, normalization, 
and scatter will require galaxy samples with uniform, well- 
defined selection criteria. 

As I discuss in §[5]below, the normalization of the ri/2— R200 
relatio n is quantita tively consistent with the expectations of 
the Mo et al. JT998) model if we assume that angular momen- 
tum of the disk is set near the end of fast mass accretion epoch 
at z ~ 2 and take int o account pseudo-evolution of halo mass 
( |Diemer et al.|2012| >. 

It is noteworthy that dlrr galaxies have similar half-mass 
radii to the dwarf ellipticals and follow approximately the 
same ri/2 - /?20o relation. Massive late type galaxies also fol- 
low the same linear relation, although the figure indicates that 
late type galaxies of intermediate stellar mass have system- 
atically larger half-mass r adii than early type galaxies of the 
same stellar mass (e.g., Berna rdTet al.|2012p . 

The shaded band around dot-dashed line in Fig. [T] shows 
2cr a; 0.3 - 0.5 dex intrinsic scatte r esti mated for all galax- 
ies in the sample of Bernardi et al. (2012 the scatter shown 
is for all galaxies in the sample, M. Bernardi, priv. comm.). 
The orange error bars show scatter estimated for the mass lim- 
ited sample of massive SDSS galaxies presented in Szomoru 



et al. ( 2012 ). The scatter estimated for this samp l e is in good 
agreement with the scatter of the|Bernardi et al. (2012} sam- 



ple. Remarkably, the scatter is also approximately consistent 
with the scatter expected from the distribution of halo spins, 
A, in models in which galaxy size is oc AR20Q, shown by the 
dotted lines in the figure. 

4.2. Stellar and gas surface density profiles of galaxies 
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Fig. 2. — Normalized surface density profiles of stars, Z„ (/?)/£„ * (top panel) 
and neutral gas, Z sas (7?)/Z„ g (HI+H2, bott om panel) for la te type galaxies. 
The top panel includes galaxies from the Leroy et al. 1 2008) and the LITTLE 
THINGS sample of Zhang et al. (2012), while the bottom panel only includes 
the surface density profiles of gas trom |Leroy et al.|(2008) . Profiles of indi- 
vidual galaxies are shown by thin lines colored according to their log 10 M, as 
indicated in the legend. Each individual profile is normalized by the radius 
r„ = 0.015R200. where R200 is obtained using the abundance matching ansatz 
from the galaxy M, . The thick lines with errorbars show the sample average 
and rms dispersion around the mean. Both stellar gas profiles are well de- 
scribed by the exponential profile on average, but the scale length of the gas 
profile is on average a factor of = 2.7 larger than that of the stellar profile. 

In this section I show that in addition to ri/2 - R200 cor- 
relation the surface density profiles of stars and neutral gas 
approximately follow universal profiles when scaled by r„ = 
0.015/?200^ i-e. the mean normalization of the r\/2 — R200 cor- 
relation. 

Two panels in Figure [2] show the surface density profiles 
of stars and neutral gas (HI+H2, where HI is corrected for 
helium) for late type galaxies. The radius of each individual 
profile was rescaled by r„ = 0.0 15/? 200 and surface densities 
were scaled by = 0.448M/>^, where M is total stellar or 
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Fig. 3. — Normaliz ed surface density pro files of stars for massive early type 
SDSS galaxies from Szomoru et al. (2012). Profiles of individual galaxies are 
shown by thin lines colored according to their log 10 M„, as indicated in the 
legend. Each individual profile is normalized by the radius r„ = 0.015^200, 
where R200 is obtained using the abundance matching ansatz from the galaxy 
M, . The thick lines with errorbars show the sample average and rms disper- 
sion around the mean. The thick dashed line shows the average profile of late 
type galaxies from the top panel of Fig.[2]for comparison. 

gas mass of each galaxy and factor 0.448 = 1.678 2 /(2tt) as- 
sumes exponential profile (r;/2 = 1.678/?^). The figure shows 
that both the mean stellar and gas profiles are on average well 
described by the exponential profile, 2(7?) = Eo exp(— /?//?</), 
where for stars So.* * 1256M*/^2oo an d me scale length 
R d , t ~ O.OIIR200 and for gas E , g a.s ~ 199M gas /7?|, and 
^d.gas ~ 0.029/?200- The gas distribution is thus on average 
a factor of w 2.6 more extended than stellar distribution. The 
scatter around the mean profiles is rather small and is only 
« 30 - 50% at R ~ 1 - 3r„, even though galaxies shown in the 
T.*(R) figure (top panel) span over six orders of magnitude in 
Af». 

The approximate universality of the gas surfac e density pro- 
files was recently pointed out by Bigiel & Blitz] p012) . These 
authors rescaled gas profiles of the THINGS/HERACLES 
galaxies using the optical Holmberg radius, R25, and surface 
density S trans at the radius where 2h 2 = Shi- Such rescal- 
ing results in average exponential profile described by £ gas = 
2.12 trans exp(-1.65/?//?25) with comparable scatter around the 
mean profile to the rescaling described above. Comparison 
with the scaling derived above gives 'Ltiws = 95Mg as /7^ m 
and/? 25 = 0.048fi 2 ()0- Thus, results of Bigi eT&^Iit^ p012) 
can be understood if 2tr an s scales with characteristic surface 
density Xo.gas- The scaling of R25 is implied by the scaling 
R d t cc l? 2 oo because for exponential disks R25 ~ 4.5/?d,»- 
Thus, the gas surface density profiles can be scaled by sur- 
face density oc M slis /R^ 5 instead of 2 tlans . In summary, results 
presented here indicate that the reason scaling employed by 
Bigiel ~& Blitz| ( f2012| ) works is that surface densities of gas 
and stars are both exponential and their scale lengths are cor- 
related: /?d,gas ~ Rd,*- The origin of the universality of the E gas 
profiles lies in the scaling of half-mass radius of both gas and 
stars with the virial radius of parent halo. 

Figure [3] shows the stellar surface density profiles of mas- 



sive SDSS (z < 0.1 ) galaxies in the sample presented by 
Szomoru et al. (2012 1 rescaled using r„ as above. Note that 
1 plot not the actual measured profiles but the Sersic profiles 
with parameters derived from M, and R e values in Table 1 in 
that paper. The mean profile of late type galaxies from the top 
panel of Fig. 2 is shown for comparison. The figure shows 
that stellar distribution of early type galaxies also follows an 
approximately universal profile. The mean profile is very 
close to the de Vaucouleurs profile with R e = 0.0 157?2oo/ 1-34, 
where factor of 1.34 converts 3D half-mass radius to the 2D 
R e . Remarkably, the mean profiles of late type and early type 
galaxies are quite similar at R > r„ and are only significantly 
different within the half-mass radius. This implies that the 
same process shapes stellar distribution at large radii in both 
late and early type galaxies. 

5. DISCUSSION AND CONCLUSIONS 

Figures [Tp] demonstrate that characteristic size of stellar 
and gas distributions in galaxies spanning more than eight or- 
ders of magnitude in stellar mass scales approximately lin- 
early with the virial radius /?200 derived using abundance 
matching approach. Such relation is in remarkable agree- 
ment with expectations of the model of Mo et al. (1998): 

n /2 oc ar 20 q. 

The agreement is not only qualitative, but also quantitative. 
The scatter in size at fixed R200 in such model is mostly due 
to scatter in A - the spin of parent halo. The dotted lines in 
Figure[T]show the scatter of 0.5 dex, which approximately cor- 
responds to 2cr lnA ~ 1.1, where c lni » 0.55 is the rms width 
of the log-normal spin pdf (e.g., |Vitvitska et al.|2002") . The 
figure shows that there is a good general agreement between 
the scatter due to the spin distribution and observed sc atter of 
r\i2 es timated for massive galaxies in the samples of Szomoru 
et al. ( 2012) ) and Bernardi et al. (2012 1. The agreement is 
noteworthy, although more detailed comparison taking into 
account measurement errors of R e and scatter of /?200 due to 
intrinsic scatter in the M, - M200 relation and measurement 
er rors of M, i s neede d- 
Mo et al.| ([1998} model predicts r l/2 = l.618R d = 



1 . 1 87 (J dlwicdfc 1 ' l fR^Rim , where jd and nid are fractions of 
baryon angular momentum and mass budget within halo in 
the central disk, f c is a function of halo concentration, and 
fit is a function that takes into account baryonic contraction 
of halo in response to halo formation (see their § 2.3). As- 
suming jdlrrid — 1, nid = 0.05, and typical spin A = 0.045 
and halo concentration of C200 = 10 gives n/2 = 0.032/?200, 
i.e. normalization about a factor of two higher than inferred in 
this study. This, of course, could s imply be due to jd/md < !■ 
However, as noted by Mo et al. (1998), the predicted rela- 
tion is applicable at the epoch of disk formation. Indeed, 
after the galaxy size is set, /?20o may increase significantly 
simply due to pseudo-evolution of halo mass defined with re - 
spect to the evolving reference density ( [Diemer et al.|2012| . 
If I assume that galaxy size was set at a characteristic epoch 
of z ~ 2 and that halo mass of most galaxies increases by 
a factor of 2.5 due primarily to pseudo-evolution between 
z — 2 and z — 0, as suggested by analysis of cosmological 
simulations ( |Diemer et al.|2012) , the halo radius that set the 
galaxy size is R2oo(z = 2) « 0.37^200(2 = 0). Using again 
A = 0.045 and C200 = 4 typical for z — 2 halos, the prediction 
is ri/2 ~ 0.016(jd/m<i)R2oo(z = 0) in agreement with the re- 
lation derived for observed galaxies. This conjecture can, in 
principle, be tested via analysis similar to the one presented 
in this paper but done at z = 1-2. 
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These estimates demonstrate that empirically derived ri/2 - 
^200 relation is in very good qu antitative agreement with pre- 
dictions of the |Mo et al7| ([1998} disk formation model. Re- 
markably, the prediction works not only for late type disks, 
but also for early type galaxies. This means that angular mo- 
mentum plays a critical role in setting the sizes of galaxies of 
all morphological types. This fact reveals yet another remark- 
able regularity in properties of observed galaxies and provides 
a critical test for models of galaxy formation. 

The derived relation may also provide a useful way to es- 
timate galaxy sizes in simulations when only halo informa- 
tion is available. Conversely, it can be used to derive halo 
extent and mass using the observed R e . As shown recently by 
Szomoru et al. (2012 1, half-light radius of galaxies is offset 
only by « 25% from half-mass radius r\ji regardless of galaxy 



stellar mass, morphology, and redshift. The relation derived 
in this study can thus be used to estimate /?200 of galaxy halos 
with w 50% accuracy and virial mass M200 to within a factor 
of about four from the measurement of half-light radius alone 
without resorting to estimate of stellar mass. For these reasons 
it would be interesting to calibrate the relation and its scatter 
using larger, well-defined samples at a variety of redshifts. 

I am grateful to Mariangela Bernardi for useful discus- 
sions and communicating scatter results for her sample and to 
Nick Gnedin and Andrew Hearin for useful comments on the 
manuscript. This work was supported in part by the Kavli In- 
stitute for Cosmological Physics at the University of Chicago 
through grants NSF PHY-0551 142 and PHY-1 125897 and an 
endowment from the Kavli Foundation and its founder Fred 
Kavli. 
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